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£SJ ■ We present a method of finding approximate analytical solutions for the spectra and eigenvectors 

of collective modes in a two-dimensional system of interacting bosons subjected to a linear external 
potential or the potential of a special form u(x, y) = fx— u cosh x/l, where fi is the chemical potential. 
The eigenvalue problem is solved analytically for an artificial model allowing the unbounded density 
of the particles. The spectra of collective modes are calculated numerically for the stripe, the rare 
density valley and the edge geometry and compared with the analytical results. It is shown that 
the energies of the modes localized at the rare density region and at the edge are well approximated 
by the analytical expressions. We discuss Bose-Einstein condensation (BEC) in the systems under 
investigations at T / and find that in case of a finite number of the particles the regime of BEC 
© ' can be realized, whereas the condensate disappears in the thermodynamic limit. 
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I. INTRODUCTION 
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The observation of BEC in alkali metal vapors confined in a trap has opened the opportunities to investigate 
experimentally this macroscopic quantum phenomenon. That outstanding discovery has stimulated theoretical studies 
of BEC in inhomogeneous Bose gases subjected to external fields (see pj| for review). One of the problems under 
, extensive investigation is the influence of an external potential on properties of low-dimensional Bose systems. It is 
' known that an external potential may cause BEC in 2D ideal Bose gases at nonzero temperatures ||-[^ . This problem 
was discussed in detail by Mullin || . Mullin has presented rigorous arguments why the Hohenberg theorem || ruling 
out BEC in 2D does not apply to the systems considered in [pH7| . 

Recently, the possibility of BEC in interacting 2D Bose gases has been investigated. Petrov, Holzmann and Shlyap- 
nikov jl0| considered a 2D Bose system confined in a harmonic oscillator potential. They have concluded that two 
regimes can be realized: a true condensate and a quasicondensate with fluctuating phase. In the mean field approxima- 
tion the properties of such a system have been studied by Bayindir and Tanatar pT[ . They demonstrate the similarity 
in thermodynamic behavior of ideal 2D Bose systems in a harmonic trap and those with weak interactions and a finite 
number of the particles. On the other hand, Mullin ||,[l2| argues the absence of BEC in trapped interacting Bose 
I ■ gases in 2D in the thermodynamic limit. 

For understanding the macroscopic quantum properties of interacting Bose systems the key point is the behavior 
of elementary excitations. In two dimensions thermally excited fluctuations destruct the condensate. Therefore, 
the information about stability of BEC can be gained from the study of the spectrum and structure of low-energy 
• • . excitations. 

The theory for elementary excitations of interacting Bose systems has been derived by Bogoliubov [^3|. After the 
experiments |l|-[| this theory has received a new interest. Stringari |l4|] has obtained an analytical solution for the 
collective mode spectrum of a 3D Bose gas in an isotropic harmonic trap in the Tomas-Fermi regime. The same 
problem for the anisotropic trap has been resolved in Ref. |L5]]. Numerical solutions of the Bogoliubov equations 
for the trap geometry have been obtained in Ref. Jpoj-pof. The exact result for the energies of the harmonic trap 
breathing modes in 2D has been found by Pitaevskii and Rosch pi] ). The collective mode spectrum in a spatially 
periodic potential has been investigated in Ref. p2| . 

Further understanding of the properties of inhomogeneous 2D Bose systems requires the study of the low-energy 
excitations in the potentials different from the mentioned above. It is important to investigate the potentials, for 
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which an analytical solution of the eigenvalue problem can be found. In this paper we analyze two such cases. The 
first case corresponds to the potential u(x, y) = fi — ucosh 2 x/l, where p, is the chemical potential. Since the external 
field pushes the particles away from the region x = 0, and a rare density valley is formed at x = 0, we call such a 
geometry an "anti-trap" one. The second potential considered is a linear one: u(x 1 y) = p, — ax. The Bose cloud 
subjected to that potential has an edge at x = with a linear dependence of the density on the distance to the edge. 

An approximate analytical solution of the eigenvalue problem is obtained by the following method. We formally 
imply the density is unbounded at — > oo (x — > +oo for the linear potential). Then, analytical expressions for 
the energies and the eigenvectors of collective excitations can be derived. To establish the applicability of the model 
description to real objects the eigenvalue problem is also solved numerically. It is shown that the spectra of collective 
modes localized at the rare density region (in the anti-trap potential) and at the edge (in the linear potential) are well 
approximated by the analytical solution. Using the analytical solution we calculate the one particle density matrix 
and analyze the off-diagonal long range order. We conclude that thermal fluctuations do not destroy completely the 
long range order in the systems with a finite number of the particles. It means that the BEC regime is realized. It is 
not the case in the thermodynamic limit, where the density of the condensate approaches to zero at T ^ 0. 

In Sec. | a general description of elementary excitations in inhomogcncous 2D Bose systems is presented. The 
collective mode spectra in the anti-trap and the linearpotential are obtained in Sec. Ill and S ec. l IV l correspondingly. 
The one particle density matrix is calculated in Sec. |v, and the conclusions are given in Sec. |vi[ 



II. COLLECTIVE EXCITATIONS OF INHOMOGENEOUS 2D BOSE GAS 

We consider an interacting 2D Bose gas in the external potential u(r). The Hamiltonian of the system can be 
written in the following form 

r % 2 » 1 » 

H = / d 2 r{ V V + (r)VV>(r) + Mr) - M+(rU(r) + -7^+(r)^+(rWrWr)}, (1) 

J 2m 2 

where tp is the boson field operator, m, the boson mass. In Eq. ([!]) we use the effective potential of the atom-atom 
interaction of the form V(r — r') = 7<5(r — r'). The strength 7 is considered as a parameter. To describe the collective 
excitations we use an approach developed by one of the authors in Refs. [P3|-p5[. We introduce the phase 0(r) and 
density p(r) operators. They are connected with the Bose field operator by the relation 

^ = e^yfp. (2) 
The commutation relation for the new operators is found to be 

p{v)(p(r') — Lp{r')p(r) — iS(r — r'). (3) 
Using Eq. (Q) we rewrite the Hamiltonian ([!]) as 

H = J d M^{^) 2 V^ + (V^) 2 + 8 [(Vv^)(V^)v^- ^)(Vv^)]}+ Mr) - ^]/5+ ^ 2 }- (4) 
The density operator can be decomposed as 

p{r) =p (r)+/3i(r), (5) 

where po (r) is determined by the condition that the expansion of Eq. (^) into powers of p\ and does not contain 
the linear terms. It yields the equation 

h 2 1 



u{v)-p + 1Po {v)-- ==V 2 V7^W=0. (6) 

2m ^p (r) 

Eq. (|^) coincides with the Gross-Pitaevskii equation for the time independent part of the condensate wave function 
Q. Therefore, the function po(r) can be understood as the density of the particles at T = Q. At low temperatures 
the fluctuations of the density and the current are assumed to be small and one can neglect the anharmonic terms. 
The quadratic part of the Hamiltonian reads as 

H 2 = [ d 2 r{-[2^f{2^) + -£Lg(-£±)] + »^-{(Vp 1 )(V^) - (V0)(Vpi) - ^[px(V£) - (V^)pr]}}, (7) 
J 4 Jpo JpH 4m po 



2 



where 

T=-^~^=i^% (8) 

G = f + 2 7Po . (9) 

Let us rewrite the operators <p and pi in the second quantization representation in terms of creation (6+) and 
annihilation (&„) operators of the elementary excitations, which should also be bosons in the system of Bose particles. 
It is convenient to write <p and p\ in the following form p4| : 



The representation (|10|), (|ll|) should transform the Hamiltonian (|7j) to the diagonal form 

H 2 = const + ^2(E V + -)b+bu, (12) 

where E v is the energy of the excitation with the quantum number v. The Hamiltonian (j?j) is reduced to the form 
(|l2|), if the functions F v (r), 9„(r) satisfy the equations 

re„ = e v f v , (13) 

G.F„ = E V Q U , (14) 

and are normalized by the condition 

J d 2 r9„(r)F;(r) = 1. (15) 

One can check that the operators 5, ip in the form @, @ with 6 V determined by Eqs. (|l|), @), and @ 
satisfy the commutation relation (g). 

One can see that under the substitution = u u + v u , 0„ = u„ — iv Eqs. (^3|), (|l^) coincide with the Bogoliubov-de 
Gennes equations We note also a similarity between our approach and the one developed by Ho and Ma [p6[ . 

III. COLLECTIVE MODE SPECTRUM IN THE ANTI-TRAP GEOMETRY 

In this section we consider the special choice of the external potential: 

u(x, y) = fj. — u cosh 2 — . (16) 

We assume the parameters of thepotential satisfy the inequality u ^> ft 2 /2ml 2 . Then, one can neglect the fourth 
term in the left-hand side of Eq. (g), that corresponds to the Tomas- Fermi approximation. At T = the density of 
the particles is equal to 

Po{x,y) = -cosh 2 -. (17) 
7 I 

To derive the spectrum of low-energy excitations (E <C u) one can omit the kinetic energy terms T in the operator 



G. Then, Eqs. (£L3|), (14) yield the following differential equation for the function O 
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2m cosh 2 j fe„(r) = E*Q„(t). (18) 



The function F is given by the expression 

F,(r) = ^cosh- 2 ^ 6„(r). (19) 
2u I 

Taking into account that the system is uniform in the y direction the function 0„ can be presented in the form 

9 v (x,y)=B v (x)e ik », (20) 

where 9 v (x) satisfies the equation 

_ tH C osh 2 j&- (fc 2 + r 2 )6 v ] = E%. (21) 
to I ax z 

Eq. ( pl| ) coincides with the Shroedinger equation for the one particle problem in the potential v{x) = cosh~ 2 {x/l). 
Its solution given in Ref. [^7| has the form 

9 u (x) = cosh-^ W ( — ^X ), (22) 



where t] = y/l + (kl) 2 , and the function w(z) is governed by a hypergeometric differential equation 

z(l - z)w" + (?7+ 1)(1 - 2z)w' - (r) - s)(r) + s + l)w = 0. (23) 



In Eq. (23) the parameter s is determined by the equation E 2 = uIqs(s + 1), where cuq = y Ti 2 u/ml 2 . To solve Eq. fl23| ) 
one should specify the boundary conditions. We consider the system infinite in the x direction. Since the fluctuations 
of the current should remain finite, the quantity d(9(x) / sj pp (x) ) Idx cannot be singular at |x| — > oo. The solutions of 
Eq. satisfying such a condition are the polynomials. Eq. j2^ ) has the polynomial solutions at s = r\ + n, where 
n = 0, 1, 2, Therefore, the energy of the collective modes are equal to 

E nk - Lo \j[n + v/T + (kl) 2 }[n + 1 + y/l + (kl) 2 }. (24) 
The exact expression for the eigenfunction 9 n k[x) reads as 

9 nk (x) = C nk cosh"" y if (tanh | ) , (25) 
where P^" is the Jacobi polynomial. The normalization factor C nk is calculated to be 

_ / MQ n!(2n + 2 V + l)r(n + 277 + 1) \ 1/2 



In Eq. (26) T is the Gamma function, L y , the size of the system in the y direction. 

The answers obtained have two uncommon features. First, the spectrum does not contain an acoustic branch. 
Second, all the modes are localized near x = 0. Such a peculiar behavior is connected with that the model contains 
an unphysical assumption. We fix the parameter I, while set an infinite size of the system in x direction. Thus, we 
allow the density of particles be unbounded at |x| — > 00. The principal question has to be clarified is whether the 
results of the model (further we call it the artificial model) can be applied to real physical objects. 

Let us consider two examples. The first one is a 2D Bose cloud of a stripe shape. We assume that there are two 
rigid walls at x — ±L X) and the potential well obeys the internal structure given at |x| < L x by Eq. ([l(]). Since the 
particles cannot flow though the rigid walls, the gradient of the phase normal to the wall should be zero at |x| = ±L X . 
Therefore, the boundary conditions for the function 9(x) are 



9(x) 



The general solution of Eq. ( pl| ) can be written in the form 



= 0. (27) 

x=±L m 
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6{x) = cosh j [A lVl {x) + A 2 tp 2 (x)}, (28) 

where 

/ n , « 1 x 1 — tanh f N 
^(2;) = cosh " 1 - aftfa - s, t? + s + 1,t) + 1, L ), (29) 

t t 1 — tsmli — 
^(i) = cosh-"- 1 T (1- tanh -)-" ^(-5,5+1,1-77, ~ 2 L )- (30) 

In Eqs. (p9[), ( p0| ) 2-F1 (a;, y, 2, u) is the hypergeometric function. We specify 1 — 77 ^ — 1, —2, .... (If 1 — 77 is negative 
integer, another form of the general solution should be used.) Using Eqs. (p7|), ( p8[ ) we get the dispersion equation 

^(i x )^(-L a ) - pi(-LxM(-&») = 0. (31) 

Its solutions can be obtained numerically. It is instructive to analyze the case of L x 3> I (at L x rs / the density 
variation across the stripe is small). The solution of Eq. ( [jl| ) at L x = 31 is shown in Fig. |l|(a). For comparison, the 
spectrum given by the artificial model (Eq. (HJ)) is also presented in this figure. One can see that in the stripe system 
the lowest collective mode is an acoustic one at small k. Its velocity is equal to c = y/jp/m, where p is the average 
density of the particles. There is a discrepancy in behavior of the two families of the curves in Fig. |l](a) only at 

E > hck. (32) 

If the inequality opposite to ( |32| ) is satisfied, one can neglect the influence of the boundaries. In this case the collective 
mode spectrum is well approximated by Eq. 

To consider another example we specify the case of a Bose cloud infinite in the x direction with the density of the 
particles given by Eq. ([l7]) at |x| < L x and the uniform density at |z| > L x . There are localized and extended collective 
modes in such a system. The energies of the extended modes satisfy the inequality (p3). The localized mode solution 



has the form 6{x) = A3 exp(Kx) at x < —L x , 9(x) = A4 exp(— kx) at x > L x , where k = \J k 2 — E 2 /h 2 c 2 . At |a;| < L x 

it is given by Eq. (p8|). The boundary conditions are reduced to the continuity of 9{x) and d(9(x) / \J po{x)) jdx at 
x = ±L X . The dispersion equation has the form 

[ip[(-L x ) - Kipi(-L x )][ip' 2 (L x ) + Kip 2 (L x )} - [(p' x {L x ) + nipi(L x )}[ip' 2 (-L x ) - Kip 2 (-L x )} = 0. (33) 

The solutions of Eq. ( |33"| ) at L x = 31 are shown in Fig. |l](b). This figure illustrates that the analytical expression ( |2^ ) 
yields a satisfactory approximation for the spectra of localized modes. From the numerical analysis we find that the 
coincidence is better for larger ratio L x /l. 

Thus, our artificial model, while is not free from some artifacts, seems quite useful for obtaining the analytical 
expressions for the spectra and eigenfunctions of collective modes for real physical objects. 



IV. COLLECTIVE MODE SPECTRUM IN THE LINEAR POTENTIAL 



In this section we consider the case of a linear external potential within the approach presented in Sec. III. This 
problem was studied preliminary in Ref. p4[ . 
The external potential has the form 

u(x, y) = p, — ax. (34) 

Let us consider the low-energy excitations in the Tomas-Fermi regime. In this regime the density of the particles is 
equal to Pq(x, y) — ax/j. 

The equation for 9(x) is reduced to 

Substituting 



%{x) = ^/2\k\xe-\ k \ x w{2\k\x) (36) 
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into Eq. (pq), we obtain the following equation for the function w(z): 

zw" + (1 - z)w' -0w = Q (37) 

with 

At the beginning, we consider the artificial model. We assume the system is semi-infinite and the density is 
unbounded at x — > oo. We are interested in the solutions finite at x = and x — ► oo. Such solutions exist if (3 = — n 
(n = 0, 1,2, . . .). Using Eq. (|38|) we obtain the following analytical expression for the energies: 



|fc|(2n+l). (39) 

m 



The function 9 u {x) reads as 



9nk(x) = C nk ^/2\k\x~e-W x L n (2\k\x), (40) 
where L n is the Laguerre polynomial. The normalization factor is equal to 



C nk = J-—, (41) 

As in the case of the anti-trap potential the spectrum of the excitations does not contain an acoustic branch. It is 
connected with the unbounded density at x — > oo implied in the model. 

Let us switch to the real systems and consider the stripe with a linear increase of the density across the stripe and 
a rigid wall at x = L x . The solution of Eq. (|35|) finite at x = has the form 

9{x) = A/x"e~ |fc| * iFi(/3,l,2|fc|x), (42) 

where \Fx{x,y, z) is the Kummer confluent hypergeometric function. Using the boundary condition Eq. ( p7| ) at 
x = L x , we obtain the following dispersion equation 

xF x {p,l,2\k\L x )-2P 1 F 1 (l + 0,2,2\k\L x )=O. (43) 

Its solutions are shown in Fig. |2|(a). The spectrum ( |39| ) is also presented in Fig. |^. In this figure the energy units 

eo = \J h 2 a/mL x are used. One can see that at k > E/hc the spectrum is well approximated by the analytical 
expression (|39|), while at smaller k the difference is essential. In particular, the lowest mode transforms to the 
acoustic one. 

To complete the picture we consider a semi-infinity system with a linear increase of the density at < x < L x and 
the uniform density at x > L x . We are interested in the edge collective excitations of such a system. The dispersion 
equation for the edge modes is found to be 

(k - |fc|) iF!(J3, 1, 2\k\L x ) + 2\k\P 1 F 1 (1 + 0, 2, 2\k\L x ) = 0. (44) 

The solution of Eq. (Q) is shown in Fig. ||(b). One can see that the spectrum of the edge excitation is very close to 
the dependence given by the analytical expression (|39|). 

To conclude this section we discuss the conditions, when the low-energy description is valid. The Tomas-Fermi 
approximation is correct at 

ft 2 

lPo(x) » J ' ( 45 ) 

when the fourth term in the left-hand side of Eq. (|^) can be neglected. In two dimensions the parameter 7 is connected 
with the radius of the interparticle interaction a by the relation [p8[ 

4-rrh 2 1 , s 

7= w 2V 46 



G 



Since Eq. ( [46| ) contains only the factor m(poa 2 ) 1 as a small parameter, the condition ( f45f ) is fulfilled at x > xq, 
where xq is of order of the local average distance between the particles. 

To omit the kinetic energy terms in the operator G we should consider the excitations with the energies 

E < ~fPo(x ) = ax Q . (47) 

Since we neglect the solution of Eq. ([57]) singular at kx — » 0, the wave vector k should satisfy the inequality 

kx < 1. (48) 

Using Eqs. (p9|), (p7[), ( fl8| ) we find the conditions on the wave vectors and the energies 

/ am \ , „ 

fc « , (49) 



E <<; Em = \ ' (50) 

Eqs. J49]), ( |50|) set the validity of the low-energy approximation. At finite L a the energies have the scale of eo- 
Therefore, the description requires eo "C E m , that yields the condition on the size of the system in the x direction 

L x > ■ (51) 

V am I 



V. OFF-DIAGONAL LONG RANGE ORDER 

It is of interest to consider the possibility of BEC in 2D systems subjected to the anti-trap or linear potential. Let 
us analyze the behavior of the one-particle density matrix 

W{r,r') = (*+(r)*(r')). (52) 

If the function W(r, r') remains finite at |r — r'| — > 00, the off-diagonal long range order is realized. 

Thermal fluctuations destruct the macroscopic quantum coherence. The principal mechanism of the destruction is 
connected with the phase fluctuations. They are described by the operator 0. Neglecting a small contribution of the 
density fluctuations we obtain the following expression for the function W(r, r'): 

-$(r, r'). 



W(r, r') « Vpo(r) Po (r0exp( ^ ; ), (53) 
where 

$(r,r')H[^(r)-^(r')] 2 ). (54) 



To be more concrete, we analyze the dependence of $(r, r') on y, y' at x = x' . Taking into account Eq. (10) we get 
$(x,y;x,y') = — - V / dk\9 nk (x)\ 2 [n B (E nk ) + - cosk(y - y% (55) 

where n B (E) = [exp(E/T) - 

Let us consider the case of the anti-trap potential. Since the maximum of the phase fluctuations is located at x = 



(see. Eq. (25)), the phase coherence destruction begins at x = 0. Therefore, if the long range order survives at x = 0, 
it takes place in the whole system. 

Using the analytical expressions (|||), ( |2^ ) for and n k(x) we obtain the following asymptotic result: 

T _l 

27rT (0) n &(0) 
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where T (x) = h 2 p (x)/2m, £t(x) — y^ 2 lPo(x)/mT 2 . In Eq. (56) the inequalities I 3> £t(0), \y — y'\ ^S> I are 
implied. The term $ zp describes a contribution of the zero-point fluctuations. This term determines the density of 
the condensate no(x) at T = 0. The temperature dependent part in the right-hand side of Eq. (pq) does not depend 
on \y — y'\. It means that the condensate survives at nonzero temperatures. According to Eq. (pq) the condensate 
density at T 7^ is connected with no(0) by the expression 

n T {0) = n Q (0)^— — J . (57) 

The result (|57j ) is obtained within the artificial model. Nevertheless, one can apply it to the real objects, if in Eq. 
(|S5|) the contribution of the modes with the energies ( |32| ) is negligible. For example, such a situation is realized for a 
Bose cloud of the shape of a square L x L with L/l 3> 1. Then, the answer ( |57| ) yields a satisfactory approximation 
for the density of the condensate. For the square geometry the total number of the particles N is connected with the 
parameter of the model by the relation 

N u ,„ I . . L. ,„ s 
Z , = ^(l + Z smh y ). (58) 



It is natural to define the thermodynamic limit as N — * 00, L — > 00 at N/L 2 = const. From Eq. (jsS) we also obtain 
L/l — const. Therefore, in the thermodynamic limit the parameter I depends on the total number of the particles 
(I ~ TV 1 / 2 ). From Eq. ( |57j ) we find that the density of the condensate decreases by the law nr(0) ~ N~ b (b > 0) 
under increase of N. It indicates the absence of true BEC in the thermodynamic limit. At finite N, as it follows from 
Eq. (pTl), the BEC regime can be realized. 

Similar calculations for the linear potential yield 

n T (x) = n (x)y - J . (59) 

In Eq. ( [sit ) the inequalities x ^> £t(x), \y — y'\ S> x are implied. For a Bose cloud of the shape of a rectangle L x x L y 
the thermodynamic limit corresponds to N — > 00, L x — > 00, and L y — > 00 at N / (L x L y ) = const and L x / L y = const. 
To analyze the density of the condensate at the same local density of the particles we should also fix the ratio x/L x 
in Eq. (B9h . Then again tit{x) approaches to zero at N — * 00 and there is no BEC in the thermodynamic limit. 



VI. CONCLUSIONS 



We have considered the low-energy collective excitations of the interacting inhomogeneous Bose gas in two- 
dimensions. We specify the two types of the external potential: the anti-trap one, in which the rare density valley is 
formed, and the linear potential, that models the simplest case of the edge geometry. We show that for the configura- 
tions considered an approximate analytical solution for the energies and the eigenvectors of the collective modes can 
be obtained. To derive the analytical solution we consider the artificial model, in which an unbounded density of the 
particles is allowed. Physically, such an approximation means that we suppress the fluctuations in the high-density 
regions. Our method yields the modes localized at the rare density valley or at the edge. Their spectra does not 
contain an acoustic branch. The modes lost in our approach can be obtained from the numerical solution of the 
eigenvalue problem. In this case we should specify the high-density profile. We present the results of the numerical 
analysis and compare them with the analytical solutions. We find that the spectra obtained numerically are well 
approximated by the analytical expressions in case of E < Tik\J jp/m, where E is the energy of the excitation, k, 
the wave vector along the valley or the edge, p, the total average density. It is important to note that in this case 
the spectrum of the collective modes is determined completely by the parameters of the external potential and does 
not depend on the interaction between bosons. Collective modes of the trapped Bose gas in the Tomas-Fermi regime 
demonstrate the same behavior juj . 

The possibility of BEC in the systems under consideration has been investigated. The high density regions affect 
significantly on the character of the thermal fluctuations in the low density areas. If the total number of the particles 
is finite, the BEC regime is realized at nonzero temperatures. This conclusion is in agreement with the results of Ref. 
fiof , where the BEC regime of a 2D Bose gas in the trapping potential has been described. The proper formulation 
of the thermodynamic limit for the systems considered shows that in this limit the true BEC does not occur, as it 
should follow from the general theorems ||^9 30 1 . 
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FIG. 1. Collective mode spectrum in the anti-trap poten- FIG. 2. Collective mode spectrum in the linear poten- 
tial for the stripe (a) and the rare density valley geometry tial for the stripe (a) and the semi-infinite system (b). Solid 
(b). Solid curves show the numerical results; dashed curves, curves show the numerical results; dashed curves, the result 
the result of Eq. of Eq. @. 
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